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Abstract 

We give a not trivial upper bound on the velocity of disturbances 
in an infinitely extended anharmonic system at thermal equilibrium. 
The proof is achieved by combining a control on the non equilibrium 
dynamics with an explicit use of the state invariance with respect to 
the time evolution. 
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1 Introduction 

In the present paper we investigate the long time behavior of an infinitely 
extended anharmonic system, that represents, of course, a schematic model 
of a crystal. There are several papers devoted to the time evolution of sys- 
tems containing an infinite number of components, either particles moving 
in a continuum or lattice systems, see Refs. [2,5,8-15,17,20,21,23], but 
few of them study the asymptotic (in time) behavior of the system. The 
reason lies in the difficulty to obtain dynamical estimates that remain good 
for very long times. Some results in this direction have been recently ob- 
tained [3,4,6,7], but they are not related to the topic of the present paper. 

Such dynamical estimates can be greatly improved if the physical sys- 
tem is assumed at thermal equilibrium. Actually, it is not necessary to 
consider exactly a Gibbs state: the present results apply to any reasonable 
time invariant (i.e. stationary) state. This assumption allows to prove the 
locality in space of the motion and hence its existence, see Refs. [1, 16, 19, 
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24,25,27-29]. Moreover, in this case it is known (but perhaps not suffi- 
ciently underlined) that the magnitude of a component of the system may 
increase only very slowly in time. This fact allows us to obtain not trivial 
results for the following problem. 

Consider a lattice system in dimension d, with anharmonic oscillators 
at each site. At time zero we perturb the oscillator located on the site i and 
study the influence at time t on that located on the site j. For harmonic 
oscillators or bounded rotators it has been proved in Ref. [22] that this 
influence becomes negligible when \i ~ j\ > ct, for some c which thus gives 
a bound of the velocity of propagation of the perturbation. 

The case of anharmonic oscillators is much more involved. The gen- 
eral case seems too difficult, while for systems in thermal equilibrium a 
reasonable estimate has been obtained in Ref. [22]: the influence of the 
perturbation is exponentially small in time for distances |« — j| > t^^^. 

At this point we must precise what we mean by perturbation. Actually, 
we do not modify the equilibrium state, rather we analyze observables 
which give an estimate of the size of the correlation. Following Ref. [22], 
we use the Poisson brackets to this purpose, by analogy with a similar 
problem in the quantum case previously studied in Ref. [18]. More precisely, 
we analyze the Poisson brackets {fi,gj o ^t} where fi,gj are observables 
localized on the site i,j respectively, and $t denotes the time evolution. We 
prove that, with probability one with respect to any reasonable stationary 
state, {fi, Qj o $j} is exponentially small in time whenever \i — j\ > t log" t 
for a suitable a > 0. This gives a non trivial bound on the velocity of the 
perturbation. 

The problems and the results are rigorously posed in Section [2] and 
proved in Section [31 

2 Notation and statement of the result 

At each point i of the d-dimensional lattice If^ there is an oscillator with 
coordinate G K and momentum pi £ R. The state of the system is 
thus determined by the infinite sequence x = {xi}igz<i = {{Qi:Pi)}i£Z'' of 
positions and momenta of the oscillators. We shall denote by X the set of 
all such possible states, equipped with the product topology. 

The time evolution t i— > x{t) = {('?i(i),?'i(0)}iGZ<^ is defined by the 
solutions of the following infinite set of coupled differential equations. 



where the force Fi{x) induced by the configuration x = {{<liTPi)}i£Z'^ on 





(2.1) 
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the i-th oscillator is given by 

F,{x)^-U'iq,)-K (ft-*). (2-2) 

j:|j-j| = l 

with > and [/(•) a non negative polynomial of degree 4 with strictly 
positive leading coefficient. Here |i — j | is the distance between the points 
i = {ii,...,id) and J = ( ji , • . • , jd) defined by 

d 

In order to consider configurations that are typical for any reasonable 
thermodynamic equilibrium state, we allow initial data with logarithmic 
divergences in the energy. More precisely, for v ^JJ^ and fc G N, let 

w,,^{x)= J2 + + E ffe-'?,)^ (2.3) 

where Ai,^k denotes the cube of center i' and side 2fc + 1. By defining 

g(x)-sup sup tI^, (2.4) 

we denote by Xq the following subset of X, 

Xo = {x e X : Q{x) < oo}. 

Let now lo be any Borel probability measure on X that satisfies the 
following superstability estimate: there exists a positive constant Coj such 
that, for any A small enough, 

u;(e^^-'=) < e^-(2fc+i)'' Vz. e Z'^ Vfc e N. (2.5) 

It is proved in the Appendix that this implies, for any A small enough, 

lim e^^w(g > AT) = 0, (2.6) 

which in particular yields lu{Xo) — 1. 

The following theorem gives existence and uniqueness of the solution to 
Eq. (|2.ip for initial data in the set Xg. 
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Theorem 2.1 There exists a one-parameter group of transformations $t : 
Xq ^ Xq, t € R, such that t — > ^t{x) is the unique global solution to 
Eq. (|2.ip with initial condition ^q{x) — x. Moreover <&t(a;)j = Xj{t) = 
{qj(t),pj{t)) are differentiable functions with respect to Xi = (qi^pi) for 
any j, i e Z'^ and i e R. Finally, there exists a constant Cq such that 

Q{Mx))^Co{Q{x)\og[e + Q{x)]+t^y (2.7) 

The proof of Theorem 12.11 is given in the Appendix. We remark that 
a global existence and uniqueness theorem for a larger class of initial data 
can be found in Ref. [17]. Unfortunately, the proof given there does not 
guarantee that the set Xq is invariant under the dynamics. We instead 
adapt to the present context the technique developed by Dobrushin and 
Fritz in the case of particles in the continuum [11, 13]. 

Let U be the algebra of all local observables. Thus f : X ^ R is 
an element of U iff f{x) = /a^xa) for some bounded set A C Z'' and 
some differentiable function /a, depending on the finite set of real variables 
XA — {xi}i£A, which is bounded with its derivatives. The action of the 
time evolution on the local observables is still denoted by : if f & U the 
function is defined by setting ^tf{x) = /(<J>t(a;)) for any x G X^. We 
finally observe that the Poisson brackets 

are well defined for any f,g^U. 

A time invariant state w is a probability measure on X for which (j2.5p 
holds and such that uj{^tf) = '^(/) for any f € U and t G R. The typical 
example of such a state is given by the infinite Gibbs measure obtained as 
the thermodynamic limit with free boundary conditions [26]. 

Given a pair of differentiable function /, g : ^ R, which are bounded 
with their derivatives, we denote by fi, resp. gj, the local observable defined 
by setting fi{x) — f{xi), resp. gj{x) — g{xj). Our main result concerns 
the asymptotic behavior of the Poisson brackets {fi, ^tgj}{x) in the limit 
when simultaneously |i — j| — > oo and < — > oo. 

Theorem 2.2 Let uj be any time invariant state satisfying (|2.5p . Then for 
each /, g as above, a > 1/2, and 6 > we have 

lim sup e'"{f„<i>tgjKx)^0, (2.8) 

*^°°i:|i-j|>tlog°t 

almost surely with respect to the probability measure oj. 
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Remark. We consider the case of one dimensional oscillators with quartic 
one-body potential and nearest-neighbor harmonic interaction only for the 
sake of simplicity. Indeed, the same strategy applies to the general case of 
TV-dimensional oscillators, say qi G M^, such that the force Fi{x) induced 
by the configuration x = on. the i-th oscillator is now given 

by 

j:|j-i|^r 

Here r is a positive parameter, while U and V are smooth functions such 
that U{^) — Pi{\S,\) and F(^) — P2(|^|), where Pi and P2 are non negative 
polynomials with maximum degree respectively 27 and ^ 7. The estimate 
l|2.8p is now valid for any a > (7 — l)/7. 

We conclude the section with a notation warning: in the sequel, if not 
further specified, we shall denote by C a generic positive constant whose 
numerical value may change from line to line and it may possibly depend 
only on the coupling constant K and the one body interaction U{-). 

3 Proof of Theorem [23] 

In this section we prove Theorem l2.2l We remark that the time invariance 
and the assumption (|2.5[) are the only hypothesis about the state uj used 
in the proof. Hence, even if oj is not a spatially homogeneous state, we 
can assume i = without loss of generality (otherwise stated, the rate of 
convergence in (|2.8p turns out to be estimated uniformly with respect to 
the location of the site i). 

Let us briefly outline the strategy of the proof. In the case of anhar- 
monic systems, the unboundedness of the Lipschitz constant of the force 
is the source of troubles in obtaining not trivial dynamical estimates. To 
overcome this problem, we consider a "good set" of initial data, the set B 
defined below, where the dynamics is under control for all sufficiently large 
times fc G N. By (|2.6p and the time invariance of cj, it is readily seen that 
uj{B) — 1. On the other hand, using rough a priori bounds on the dynamics 
for short times, an obvious interpolation shows that in the set B the time 
evolution $4 is under control for any t ^ 0. Such a control turns out to 
be good enough to solve and bound, by iteration, the (linear) variational 
equation for the disturbance. We now proceed with the proof. 

Given a > 1/2 as in the statement of Theorem I2.2[ we choose 5 G 
(1, Aa — 1) and define 

CXD 00 

B=\J f]Bk, Bk^{xeX : QiM^)) ^ log' k}. 

n—1 k—n 
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Setting Bl = X\Bk, since lu is time invariant, we have 

uj{bI) =cj(Qo$fe >log'^fc) ^u;{Q >\og^ k). 

It follows, by dMl), that if A is small enough then e^^°i^' ''u;{bI) ^ as 
k oo. In particular, since 6 > 1, J^k ^i^k) < whence 



u{B) = 1 



(3.1) 



by the Borel-CantcUi lemma. 
We next estimate: 



dqo dpo dpo dqo 



dfo f dgj 
dqo V dqj 



lit) dgj 



dqjjt) ^ dgj 
dpo \ dqj ^^(^) dqo dpj 

^ 4||V/||oo||Vg||oo||A,(t,x)||. 



dpo 
dp J it) 
<i>t(x) dqo 



(x) 
(3.2) 



Here {qj (t) , pj (t)) denote the j-th coordinates of <I>f(x) and ||Aj(t, a;)|| is 
the uniform norm of the 2x2 Jacobian matrix Aj(t, x) given by 



Ajit,x) 



d^tix)^ 
dxo 



I dq,{t) dqjjt) \^ 

dqo dpo 

dpjjt) dpjjt) 

\ dqo dpo / 



(3.3) 



Hence, by (|3.ip and (|3.2p . Theorem 12.21 follows once we show that, for any 
6> 0, 

lim sup e^*\\Aj{t,x)\\=0 'ixeB. (3.4) 

*^~i:|il>tlog°* 



Lemma 3.1 There exists a positive constant Ci > such that, for any 
X e Xo, j G Z"*, and t > 0, 



\Ajit,x)\\ ^(l+t) 



H{t,x) log^/\e + n)Y 



where 



H{t,x)^Cit^ sup Qi^six)) 



1/2 



(3.5) 



(3.6) 
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Proof. By (|2.ip . the trajectory {qj {t) , pj (t)) = ^t{x)j satisfies the equation 

Mt))'[ P3 ) Jo \ Pji^M) 

from which we get, recalling (|2.2[) . 



is,x), (3.7) 



with 



= -[[/" (g,(s)) + 2dif] 5,-;, + J2 ^^'h. (3.8) 

f:|£-j| = l 

The integral equation (|3.7p can be solved by iteration, getting 



1 t 



Aj{t,x) = (q 1 ) f^oj + 2^ I dti I dt2--- I dtnGjih, . . . ,tn) 



(t - ii)(ti - ia) • • • - tn) 0\fl t„ 

[tl ~ t2) ■ ■ ■ {tn-l - tn) OMO 1 



(3.9) 



where 



5o,fc,(ii)---Bfe„j(t„) if lil 

Gj{ti, . . . , tn) ~ < fci,...,fc„ 

otherwise. 



From the definition (|3.8p . by using (|2.3|) . (|2.4p . and recalling t/(-) is as- 
sumed a non negative polynomial of degree 4 with strictly positive leading 
coefficient, it follows that 



\Bk.h{s)\ ^ CVQ($s(x)) log(e + |fc|) VA:,/ieZ'* VseR. (3.10) 

We now observe that the sum in (|3.9p involves only sites kg such that 
^ n. On the other hand, the number of n-step walks, starting from a 
given site and with the possible presence of some permanences, is bounded 
by (2fi-|- 1)". Then, by ^ and (IHT^ . 

iiA,(t,.)|i ^ ii+t) y: . 

By the Stirling formula the bound (|3.5p now follows for a suitable choice 
of Ci in 
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We can now prove (|3.4p . Let x G B, hy the definition of B there exists 
a positive integer uq = no{x) such that Q{^k{x)) ^ log'' k for any k ^ no. 
By (I^TT)) . for any t > no, 

sup Q($5(a;)) < sup Q{<^s{x)) + max sup ($fc(a;))) 

< Co{Qix)\og[e + Q{x)]+4} 

+ max Co{g($fc(x))log[e + g($fe(a:))] + l}, 

fc=no, ...,[*] 

whence, for Iq = io(a;) sufficiently large, 

sup g($s(a;)) Clog'^Hoglogt \ft ^ io- 

By p.Sp and p.6p we thus obtain, for any t large enough, 



\A,{t,x)\\ + ^ 



Ct^ J log(e + 7i) log" t log log t 



If |j| > i log" t, recalling S < 4q; — 1, the right hand side is readily seen to 
be bounded by exp { — t log^^^ t} for t sufficiently large. The limit (|3.4p is 
thus proved. 

Appendix 

In this appendix we prove Eq. (|2.6|) and Theorem 12. II 
Proof of Eq. (j2.6p . From the definition (|2.4p we have 



ioiQ>N)^J2 J2 to{W,^k> N{2k + l)'^). 

Applying the Exponential Chebyshev Inequality and the assumption (|2.5 
on the measure u, we next bound, for some Ao > sufficiently small, 

Lo{W,^k > N{2k + if) ^ ^(e^o(H/„,.-Ar(2fc+i)'^)^ ^ ^-(NX,-c^)[2k+ir _ 
Then, for any N > {d + 2C^)/Xq, 

-NXo(2k+lY /2 

lyel."^ fc>logi/'*(e+k|) 



c.(Q>iV) ^ 5] ^ 
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which implies (|2.6p for any A < Aq 



Proof of Theorem 12. li For any initial condition x = {a^ijieZ'' ^ ^o, the 
solution to Eq. (|2.ip is constructed as the limit 



= lim <i>'r\x), VzeZ^ 

n— >oo 



(A.l) 



where the n-partial dynamics ^["'\x) is defined in the following way. For 
any n e N let A„ = Ao,„ be the cube of side 2n + 1 centered in the ori- 
gin. Then $i"^(a;) = {(g-"''(0;-P!"^(0)},gA is the solution to the Cauchy 
problem 



t(g(")(0),p(")(0))=x„ 



i e A„ 



where 



ieA„:|j-i| = l 



(A.2) 



(A.3) 



Lemma 3.2 There exists a positive constant C2 such that, for any x G Xq, 
t G K, n G N, and Q A„, 



W,^k{^\''\x)) ^C2{Q{x)[\og{e + n) + k''] +t^), 
where Wi,^k{') is defined in (|2.3p . 



(A.4) 



Proof. For notational simplicity we consider the case i > 0. By the 
equations of motion, 

where J2* denotes the sum over all the sites i G A^^k and j G A„ \ A^^s, 
such that \i — j\ = 1. Recalling the assumptions on the one body potential 
[/(•), we have 



dt 



1/2 t 



1/4 



On the other hand, Wu,k+i ^ Y^'^i W^'^k, where Y^'^, denotes the sum over 
all the sites v' such that — ^ 1 for some £ = 1, . . . ,d; then, letting 

Wk{t)^ max 

iy:Ai,.fcCA„ 
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we get 

d 



Now, by integrating the last inequality from to t and taking the maximum 
for v such that A,y,fc C A„, 

Wk{t) < Wk{0) + C f dsWk{sf/\ 



whence Wk{t) ^ (lyfe(O)i/'* + Ct)^. But from the definition (EH) we have 

that Wfe(O) < Q(a;){2[log^/''(e + n) + k] + 1^. The bound (TOl) follows 
from the previous estimates. ■ 

We now prove the existence of the limit (jA.l[) . Again we consider the 
case t > 0. We define: 

Uk{n,t) = max(5i(n,t), k ^ n, 

leAk 

d„(t) = max max{|g(")(.)-gf)(0)| + |p(")(s)-p("H0)|}. 
From (jA.2p and (jA.3p . for any z G A^, we have 

Jo 



+ / ds|F,($("+i)(a;))-F,(ci>(")(a;))| 
Jq 

^ f ds[l + 2dK+\U"i^^{s))\]uk+iin,s), 



where Cj"(s) is a point in the interval with endpoints ql"\s) and ^^"^^^(s). 
By (TO)) , for any s e [0,i], 

|C/"(C(s))| ^ ^^[1 + Qt\s)' + ql"+'\sr] ^ C{Q{x) log(e + n) + t^}'^', 

so that, setting (pn{t,x) = {Q{x) log(e + n) and taking the maximum 

for j £ Afc, we arrive at the following integral inequality, 

Uk{n,t) s^C tp„{t,xY/'^ I dsuk+i{n,s), 



which can be solved by iteration getting, for any n > k, 

j-n — k 

Ukin,t) ^ C"-'=^„(t,:E)("-'=)/2 r^d^it). 

[n — k)\ 
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By using the Stirling formula and observing that 

d„(t) < t max max{\p["\s)\ + Ct^„(i,x), 

we finally obtain 

Mn,t) ^ l^^^J^k ■ 

Choosing nfe = nk{t,x) ~ 2k + C3[l + Q{x) +t\* with C3 large enough, it is 
easily seen that Uk{n,t) ^ 2~*^"~'"') for any n ^ Uk- It follows that Uk{n,t) 
is /c-summable, which implies the existence of the limit (jA.l[) . Indeed, we 
have also an estimate on the rate of convergence: for some constant S > 0, 

- <^>["-\x),\ e-'"'' VieAfc. (A.5) 

The proof of the differentiability of $t {x)j with respect to Xi is quite 
standard and we omit the details. Note however that we explicitly solved 
the corresponding variational equation (for i — 0), see Eqs. p.7p - (|3.9p . 

Finally, we prove the estimate (|2.7p . Obviously, this also shows that 
^t{x) e for any t g R. Given 1/ G Z"^ and A; G N, we set = x) 
and estimate 

W,,k{Mx)) ^ W^-,fe(*i"*^(^)) + \W.,k{M^)) -W,,k{<i>'t"'\x))\. (A.6) 

By (Q, for any k > log^/''(e + \iy\), 

W,^k{^'i'"'\x)) C2{Qix)[\ogie + n,) + k'']+t^} 

^ Ci2k+lf{Q{x)log{e + Qix))+t^}. (A.7) 

On the other hand, 

X max _pf + - Uiqi'''\t)) 

< C{2k + l)\„,{t,xf^^ max - 

ieAi,,fe 

Ci2k + lfip^,{t,xf/*e^^''' «C C, (A.S; 
where we used ()A.5p . From (|A.6|1 . (|A.7p . and (jA.SP the bound (pjj) follows 
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